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TECHNICAL NOTE 2822 

A SPECIAL INVESTIGATION TO DEVELOP A GENERAL METHOD FOR 
THREE-DIMENSIONAL PHOTOELASTIC STRESS ANALYSIS 
By M. M. Frocht and R, Guernsey, Jr. 


SUMMARY 


The method of strain measurement after annealing is reviewed arid 
found to be unsatisfactory for the materials available in this country. 
A new, general method is described for the photoelastic determination 
of the principal stresses at any point of a general body subjected to 
arbitrary loads. The method has been applied to a sphere subjected to 
diametral compressive loads. The results show possibilities of high 
accuracy. 


INTRODUCTION 


It is known that purely photoelastic procedures cannot solve the 
general three-dimensional stress problem. The photoelastic method fur- 
nishes five independent equations, whereas the complete specification 
of the state of stress at a point requires six relations to determine 
six unknown stress components. 

In order to obtain a sixth relation it has been suggested that the 
frozen slices removed from the model be annealed and strain measurements 
be made after annealing. This suggestion has recently received a rather 
extensive treatment from Prigorovsky and Preiss in Russia (reference 1). 

A careful analysis of this suggested method shows that its successful 
application requires model materials having relatively low values of 
Poisson's ratio at the elevated temperatures used in the freezing proc- 
ess. Such materials are not available in this country. Fosterite and 
Bakelite, which are the best available materials, have Poisson's ratios 
approximately equal to l/2. It is further shown that the method of strain 
measurement after annealing bre.aks down when this ratio approaches l/2. 

In this report a new method is described which does not depend on 
Poisson's ratio and therefore can be used with models made of Fosterite 
and Bakelite. This method employs frozen stress patterns from normal ■ 
and oblique incidence. The separation of the principal stresses is 
obtained by the mmierical integration of one of the differential 
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equations of equilibrium in Cartesian coordinates rather than by strain 
measurement after annealing which involves Poisson's ratio. It will be 
shown that this .permits the determination of all six stress components 
at each point of a body. 

The report consists of three parts. . The first part comprises a 
survey and analysis of the method in three-dimensional photoelasticity 
which rests on the freezing and slicing processes and strain measure- 
ment after annealing. The second part presents the theory of the new 
method. The third part contains the application of the new method to 
the determination of stresses in a diametrically compressed sphere. 

The investigation was conducted in the Photoelastic Laboratory of 
the Mechanics Department at the Illinois Institute of Technology under 
the sponsorship and with the , financial assistance of the National Advisory 
Committee for Aeronautics. ' 

The Research Corporation provided the funds for the fellowship held 
by Mr. Roscoe Guernsey, Jr. 

Mr. David Landsberg, Assistant Research Engineer in Experimental 
Stress Analysis, assisted in all experimental phases of the work. It 
is a pleasure to acknowledge his cooperation. 

Acknowledgment is also due Mrs. Dora L. Frocht for her assistance 
in the translation of the paper by Prigorovsky and Preiss (reference 1 ). 


SYMBOLS 


^x^ cfy# 

normal stresses, psi 



’’’xy^ "^yz^ ■’’zx 

shearing stresses, psi 




stresses on an inclined plane 

, psi 


Pf q 

principal stresses, psi 



pS q* 

secondary principal stresses 

in XY -plane. 

psi 

p", q" 

secondary principal stresses 

in XZ -plane. 

psi 

Oyf 1^2 

normal stresses in fringes 



“xy> “yz^ ^zx 

shearing stresses in fringes 
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r, 0 ", f' 


0ey 

F 

F’, F" 


f 

E 

V 

P 

A 

Ac 

R 

Rc 


T T . T 

xy^ yz^ zx 


«yz 


fringe orders at normal incidence for slices parallel 
to XY -plane, XZ-plane, and YZ-plane, respectively 

fringe orders at oblique incidence for rotations about 
Y-axis and Z-axis, respectively 

angle of rotation of a slice j also angle defining an 
inclined plane 

isoclinic parameters at normal incidence for slices 
parallel to XY -plane, XZ-plane, and YZ-plane, 
respectively 

isoclinic parameter at oblique incidence for rotation 
about Y-axis 

shear fringe value of slice, psi per fringe 

shear fringe value for slices parallel to XY -plane and 
XZ-plane, respectively 

shear fringe value of material, psi per fringe per inch 
linear strains, inches per inch 

Young’s modulus, psi 
Poisson’s ratio 
load^ pounds 

area of equatorial plane of sphere, square inches 
area of surfaces of contact, square inches 
radius of sphere, inches 
radius of contact areas, inches 
normal stresses in terms of P/A 
shearing stresses in terms of P/A 
contact pressure in terms of P/A^ 
contact shearing stress in terms of P/A^ 


t 
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SURVEY AND ANALYSIS OF EXISTING METHODS 
Frozen Stresses and Oblique Incidence 


Frozen pattern ,- It is now well established that elastic stress 
systems can be fixed or frozen into models made of certain diphase 
plastics and that such models with frozen stresses can be sliced into 
thin sections without disturbing the fixed pattern (references 2^ 3 ^ 
and 4), Observations of such slices in a polari scope yield the rela- 
tive retardations as well as the isoclinic parameters at each point of 
the pattern. 

Oblique incidence ,- The use of oblique incidence of a collimated 
beam of light, as suggested by Drucker and Mindlin, adds materially to 
the infonnation obtainable photoelastically (references 5 6). The 

retardation and isoclinics at normal incidence are a function of the 
secondary principal stresses in the plane of the slice, while those at 
oblique incidence depend on the secondary principal stresses in a plane 
perpendicular to the wave normal at each point of the slice. 

The basic relation for oblique incidence with rotation about the 
Z-axis is given by the following expression 


cos^0 


( 1 ) 


. The system of notation used in this report is shown in figure 1. 
Normal stresses are positive when tensile and negative when compressive. 
The four components of shear in the XY -plane are referred to either as 
the Txy or the Ty^ shear system, and the sign of this system is 
positive when the shearing components are as shown in figure 1 (refer- 
ence T, par. 1.3). Similarly the shear system in the YZ-plane is posi- 
tive if the components are as shown in figure 1. No signs are attached 
to individual shearing stress components, their directions being deter- 
mined by inspection (reference 7^ pax. 8.2). 

By combining the data from five stress patterns of different, 
obliqueness it is possible to determine the three differences between 
the normal stress components and the three systems of shearing stresses 
at each point in the slice (for convenience the plane of the slice is 
taken as one of the coordinate plane s ) . 


^("’’xz cos 6 + Ty 2 sin 6) 


2 




(cTx - 1^2) + (oy - Ox) sin^0 + T 


xy 


sin 20 


S' 
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It can be shown that from the five quantities obtained with the 
aid of oblique incidence it is possible in turn to obtain the three 
principal shears at all points of the section. This is equivalent to 
determining Mohr's circle for a three-dimensional state of stress except 
for its origin which remains indeterminate. 

Limitations of purely photoelastic data .- Except for special cases, 
the optical data by themselves are insufficient for the determination of 
the individual principal stresses. This limitation results from the 
fact that Isotropic stress systems produce no photoelastic effects. Con- 
sequently, two states of stress differing by an arbitrary isotropic 
system produce equal photoelastic effects. 

The method employing scattered light, or the Tyndall effect, which 
was developed in this country by Weller (references 8 and 9) and inde- 
pendently by Menges (reference 10) suffers from the same limitation. 

The method of convergent light employed by Hiltscher (reference 11) 
and by Kuske (reference 12) makes it possible to determine also the 
directions of the principal stresses but not their magnitudes. 


Separation of Principal Stresses 

Free surfaces .- The limitation mentioned above does not hold at 
free boundary surfaces. A free surface is subjected to only two prin- 
cipal stresses, similar to those in plane stress systems. Tangential 
slices yield directly the difference between the principal surface 
stresses. If in addition a slice is tak:en normal to the surface and 
parallel to one principal stress, it is possible to determine the indi- 
vidual principal stresses on the surface (fig. 2), This method has been 
employed by Leven and Frocht (reference 13) to determine the principal 
stresses on the surface of Diesel engine valves. Leven (reference ll<-) 
has also applied this method to the problem of surface stresses in tor- 
sion, and Hetenyi (reference I5) has applied it to threaded connections. 
In these applications the faces of the slice were oriented to be normal 
to the direction of a collimated polarized beam. 

In special cases the combination of oblique and normal incidence' 
leads to a complete determination of the principal stresses. Using this 
combination, the stress distribution in Saint Venant torsion was deter- 
mined (reference I6). 


Planes of symmetry .- For the special case where a plane of symmetry 
exists Jessop (reference 17) has developed ein extension of the Lame- 
Maxwell equations (Filon’s graphical integration) to three-dimensional 
cases. By means of these "extended equations it is possible to determine 
the stresses along the axis of symmetry. However, the method lacks 
generality. 
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Strain measurement after annealing .- It has been suggested by Kuske 
(reference 18) that mechanical strain measurements after annealing in 
conjunction with the freezing method might be used to provide the addi- 
tional relation necessary for the determination of the principal stresses 
at a general point. If it be assumed that the differences between the 
three normal stress components at a point have been found photoelasti- 
cally from equation (l), there result: 



( 2 ) 


where the C*s represent constants. If now the slice from the model 
containing the frozen stresses is annealed, the state of stress i^ 
relieved and the thickness of the slice at each point will return 'to 
its original unstressed value. From this change in thickness, if it 
can be measured with accuracy, the strain at a point in a direction 
perpendicular to the slice can be computed. Taking this direction as 
the Z-axis, for instance, the strain would be €2. Then from Hooke’s 
la\/> 


€z = ^[5^ ~ v(ox + Oy^ 


(3a) 


which may be written as 


€z = ^ ~ ^y) + '^(®z “ (3t>) 


from which 


0 > 7 . — 


- ‘^(oz ~ t^y) ~ ^(^z ~ Qx) 
(1 - 2v) 


(ka) 


In view of equations (2) this gives the stress component Cz after 
which Ox and Oy are readily found. . The entire state of stress has 
thus been determined. 
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Limitations of Mechanical Strain Measiirement 

Poisson* s ratio equal l/2 .- The method outlined above would seem 
to solve the problem and offer a powerful method of attack. Closer 
examination discloses certain serious difficulties. For most of the 
plastics used in this country Poisson's ratio is very nearly l/2 at the 
elevated temperatures used in the freezing process. If Poisson's ratio 
is taken as l/2 then inspection of equation (Ua) shows that the denomi- 
nator vanishes and the equation is quite useless. The method evidently 
fails xinder these circiomstances. This is a manifestation of the general 
proposition from the theory of elasticity which states that when Poisson's 
ratio is l/2 a given strain field does not determine a unique stress 
field, although the converse is not true. A given stress field does 
determine a unique strain field. 

Thus, when Poisson's ratio is l/2 the general equations of Hooke's 
law may be written in the form below: 


Ox = 


~ E€x + CTi 


= 3 ^^y + 


Ov. — 


^ Ecz + ai 


(5a) 


where 


0^ + a + 0^ 

0i = — ^ 

^ 3 


For this set of equations the determinant of the coefficients vanishes 
and the matrix is of rank two. There are, therefore, infinitely many 
solutions. Furthermore, these solutions differ from each other by arbi- 
trary isotropic stress systems 0 j^. 

As already stated, in the method employing strain measurement after 
annealing the procedure is to obtain photoelastically two expressions 
giving the differences between the normal stress components and to combine 
these with one expression based on the measurement of one strain nor- 

mal to the face of the slice. This yields the following equations: 
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a 


z 


'^x “ ^zx 


N 



(5b) 


As in the case of equations (5a) the determinant of the coefficients of 
equations (5b) also vanishes and the matrix is' of rank two. There are, 
therefore, infinitely many solutions to equations (5b). Consequently, 
the method of strain measurement after annealing fails when v = l/2. 

The situation is exactly analogous to that existing in photoelasticity 
itself where the principal stresses at a point can be determined only 
up to an arbitrary isotropic system. Strain measurement when v = I/2 
thus adds nothing to the general solution of the three-dimensional 
problem. 

, Poisson’s ratio nearly 1/2 .- If Poisson's ratio is slightly less 

than 1/2, the method of strain measurement after annealing should theo- 
retically lead to a solution. However, other difficulties arise. Equa- 
tion (4a) may be written as 


Ov. — 




(1 - 2V) 


(4b) 


If V is the true value of Poisson's ratio and v + Av is its experi- 
mentally determined value, then the error Aa^ in the computed stress 
Oz* = ®z + S’ measured eg is 


tJz ” o’z = - 


Av 


1 - 2V - 2AV 


(cTx + cry) 


(6a) 


With V very close to I/2, a very small error AV in Poisson's ratio 
may lead to large errors Aa^ in the con^juted stress. For instance, 
if V = 0.48 is assumed, which is the approximate value for Fosterite, 
and Av is taken as only 0.01, then 

^CTz = Q*Q 2 ^^^ “ 0.50(Ox + CTy) 


(6b), 
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Experimental measurements of v »- Experiments with the determination 
of Poisson’s ratio for Fosterite and similar plastics indicate that it 
will be rather difficult to determine the value of V closer than 
±5 percent.- The error Auz in equation (6b) would be particularly large 
when Ox and Oy happen to be of the same sign and each is large in 
comparison with cr^. All things considered, no great accuracy can be 
expected from this method so long as Poisson’s ratio is nearly l/2. 

It must be pointed out, however, that strain measurements may serve 
a useful purpose. Assuming that, in some way or other, the normal stress 
components have been found, the strains can be calculated and compared 
with those found experimentally. Here the error in the computed strains 
due to an error in Poisson’s ratio is given by 


A^z 




X + 


(T) 


which is not large. 


Method Suggested hy Prigorovsky and Preiss 

The method outlined above for the separation of the principal 
stresses which employs oblique and normal incidence of collimated polar- 
ized light and strain measurement after annealing is not the only possi- 
ble procedure. Prigorovsky and Preiss suggest two alternative methods 
in reference 1. Their procedures combine (l) stress patterns from nor- 
mal and oblique incidence with (2) axis patterns from convergent polar- 
ized light and (3) strain measurements after annealing. The significant 
point lies in the fact that the final equations obtained from their pro- 
cedures have identically the same limitations as equations (5)* 


THEORY OF SHEAR DIFFERENCE METHOD 


General theory .- A method for determining stresses in three- 
dimensional problems is now proposed which is completely general. With 
this method the six stress coii5)onents at any point may be found. It is 
essentially an extension to three dimensions of the method, long and 
effectively used for plane problems, which is known as the shear differ- 
ence method (reference 7^ ch. 8). 
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Consider an arbitrarily loaded unsymmetrical model with the set of 
coordinate axes, as shown in figure 3 . Let a straight line AB be drawn 
through point (i) from boundary to boundary and let this line be taken 
as the X-axis. At any point along this line the first differential equa- 
tion of equilibrium, with body forces neglected, is 


— + ^ = 0 
8x by bz 


( 8 ) 


and upon integration the stress at any point (j) is given by 


((Tx)j = (ox)a - 


8 t- 


22E 


8y 


dx 


-i: 




zx 


8z 


dx 


( 9 ) 


where (oTx)a denotes the stress at point A and "the stress at 

any other point (j) on the line AB. The partial derivative — is the 

8 y 


rate of change of Tyx with respect to y and 


'zx 


bz 


is the rate of 


change of with respect to z. Thus if values of Ty^ were computed 

along a line through (i) parallel to the Y-axis and the curve xyx = f(y) 


were 


plotted, then would be the slope of the curve f(y) at 

( 4 ), ■' 


point (i). Similarly 


is the slope of the curve tzx = s(z) ' 


at point (i). As in the plane problem, these slopes may be approximated 
by computing the shearing stresses at points near (i) on opposite sides 
of the line AB and forming finite difference quotients. Thus, choosing 
nei^boring points C and D in the XY -plane and s imi larly points E and F 
in the XZ -plane. 


^ 8'*'y:A ^ ('^yx)c ~ ('''yx)D _ /A'’'yx\ 

^ 8y /i ’ Ay , ■ \ Ay /i 


(10a) 
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^^zx\ ^ (^zx^F ~ 

hz./i Az 


/AT 


ZXl 


Az 


(10b) 


Thus if the shearing stresses can be determined along four auxiliary- 
lines, parallel to and on opposite sides of AB, two lines in the XY -plane 
and two in the XZ -plane, one has all the data necessary to obtain the 
quotients’ on the right side of equations (lO) and hence good approxima- 
tions to the partial derivatives. In evaluating equations (lO), care 
must be taken to attach the proper signs to the shear systems Tyx and 
Tzx, as in figure 1. 

Substituting the above approximations for the partial derivatives 
in equation (9) and replacing the integrals by summations, the following 
equation is obtained: 


('"x)j = ('^x)a 



^yx 

Ay 


Ax 



(11a) 


The summations are evaluated graphically in the same manner as in plane 
problems. For convenience. Ay and Az may be taken equal to Ax. 
Then. equation (lla) becomes 

6 j 

(<Jx)j = (Tx)a - "S ^Tyx ~ ^ At^ (11b) 

a a 


in which At,„ and At have the mean value in each interval Ax. 

Shearing stresses in first slice .- In order to carry out this 
integration^ it is necessary to determine the magnitude and direction 
of the shearing stresses along the tvo auxiliary lines in the 

XY -plane and of r^x along the two auxiliary lines in the XZ-plane. 
The shearing stresses Tyx are obtained from a slice in the XY -plane 
containing AB in its middle surface. The stress pattern of this slice 
from normal, incidence will give the difference between the secondary 
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principal stresses in the plane of the slice at all points, and the 
corresponding isoclinics furnish their orientation. The magnitude of 
the shearing stress Ty^ at any point will then be given by 


T 


yx = |(p' - q') sin 20* 


• (12a) 


where p* and q* are the secondary principal stresses in the XY -plane 
and 0’ is the isoclinic parameter.' The directions are determined hy 
inspection as in paragraph 8.2 of reference 7- Using equation (l2a)^ 
the shearing stresses Tyx along the auxiliary lines and along AB itself 
may be found. 

Shearing stresses in second slice .- A second slice lying in the 
XZ-plane and containing line AB would furnish similar information for 
Here a practical difficulty arises since the first slice removes 
an essential part of the second slice. One of several procedures may 
be used to eliminate this difficulty. 

(1) In the general, case two identical models, identically loaded, 
may be used, one for the XY slice and one for the XZ slice. The shearing 
stresses xyx for the XY slice are calculated from equation (l2a). Simi- 
larly, the shearing stresses zx. for the XZ slice are given by 

T2 x = |(p" - q") sin 20" (l2b) 

where p" and q" are the secondary principal stresses and 0" is the • 
isoclinic parameter in the XZ slice. 

(2) In large models it may be possible to use a sub slice from the 
main, slice for determining t^x.^ After the data are obtained from the 
main slice, a small section containing the line AB is cut from it, as 
shown in figiire 4(a). The necessary data in the XZ-plane are then 
obtained from normal incidence on the sub slice, as shown. This procedure 
is feasible if the model is large so that the main slice can be made of 
sufficient thickness. 

« 

(3) In the particular case where a plane of stress symmetry exists, 
such as the XY -plane in figure 4(b), advantage can be taken of this 
symmetry. Referring to figure 4(b), let it be assiuned that the stresses 
on line AB are required. The first slice is made parallel to the 
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XY -plane and contains the line AB as shown in figure 4(c). The orthogo- 
nal slice is cut from the opposite side of the body so that it contains 
the symmetrically placed line A'B', along which the stresses are the 
same as ailong AB itself. The necesseiry shearing stresses are calculated 
as outlined in procedure (l), above. 

Normal stresses .- The starting value (cTx)a will be obtained from 
boundary conditions and boundary fringe orders. The integration may 
then be carried out and values of Ox obtained along AB. Further, from 
Mohr's circle or other considerations; 

(ax - ffy)j = (p! - q')j cos 20'j = 2F'n'j cos 20'j (l3a) 


(<Jx - = (p" - q")j cos 20". = 2F"n"j cos 20"j (I3b) 


where the F's and n's denote, respectively, fringe value of the model 
in shear and fringe order at point (j). From equations (13a) and (l3b) 


(^y)j = (®x)j - (p* - q')j cos 20 ' j 
(®z)j = (<^x)j - (p" - q")j cos 20 "^ 

All the necessary data for the evaluation of 

obtained from the slices in the XY- and XZ-planes, respectively. 

Use of oblique incidence .- At this stage five of the six stress 
components, namely Ox, Cy, Oz, Tyx> and have been found at 

all points of AB. There remains one unknown stress component Tyz> The 
shear system Tyz has no influence on the stress patterns from normal 
incidence for either one of the slices but will have an effect on the 
patterns from oblique incidence. 

In order to find the shear system Ty^, an oblique stress pattern is 
obtained from either one of the two slices. For concreteness assume that 
the slice parallel to the XY -plane is used and that it is rotated in a 
clockwise direction about the Y-axis through an arbitrary angle 9. 


( 13 c) 

(13a) 



Figure 5 shows a view of a small element as seen from the positive end of the Y-axis. From 
equation (l) the fringe order ngy at any point due to the oblique incidence is given by 
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and 


sin 200y = 


2(nyz cos 9 + nyx sin 0) 
ngy cos 0 


(15b) 


Similar equations may be written for rotation of the slice parallel to 
the XZ-plane about the Z-axis. 

It is to be noted that in general the retardation observed at any 
point depends on the direction of rotation of the slice. For one direc- 
tion of rotation the fringe order at a point will he different from its 
value for the other direction. In the particular case when the slice 
contains a principal plane, then = T^y = 0 eind equations (l4a) and 

(lUh) become identical. In such cases the direction of rotation is imma- 
terial. In dealing with general slices it is important to note carefully 
the direction of rotation relative to the wave normal and to attach the 
proper signs to all the stresses. 

Each of equations (l4a), (lUb), (15&)^ an<i (l5b) may be solved for 
the unknown, shear system nyz. It is necessary only to determine the 
fringe order and the isoclinic parameter 0Qy along the line AB. If 
the rotation is counterclockwise equation (l4b) or (l5b) is appropriate. 
Of these, equation (l5b) is much the simpler. Using equation (l5b) and 
0 = 45° there is obtained 


nyz = I n^y sin 200y - nyx (l5c) 

With this, and = 2Friy.^ are easily computed. All six com- 

ponents of stress are thus determined for the point (i), and therefore 
the principal stresses themselves are determined at the given point. 

Stresses in the plastic range .- It should be noted that the method 
described in this report is not limited to the elastic state. It is 
equally valid in the plastic state of the model, provided the stress- 
optic law for this reinge is known. This follows from the fact that the • 
only equations, in addition to the stress-optic law, used in this method 
are the equations of equilibriiim, which are independent of stress-strain 
relations. It should, however,, be noted that the results will apply to 
the model only and that they are not directly transferable to the 
prototype . 
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Effect of Poisson's ratio .- In conclusion, it should be noted that 
in transferring the results from three-dimensional photoelastic 'models 
to metal prototypes the effect of Poisson's ratio will have to be con- 
sidered. It is fortunate, as shown by the theoretical solutions obtained 
to date, that the effect of Poisson's ratio on the most significant 
stresses is small (references I9 and 20). 


APPLICATION OF SHEAR DIFFERENCE METHOD TO 
A DIAMETRICALLY COMPRESSED SPHERE 
Description of Apparatus 


The apparatus used in this investigation consisted of the following 
items : 

(1) An electric furnace with temperature controls and built-in 

loading frame with special jigs 

(2) An 8-inch photoelastic polariscope with a special immersion 

tank ■ . 

(3) An oblique incidence jig 

(4) A Babinet-Soleil compensator 

A photograph of the electric furnace is shown in figure 6. This is a 
relatively large oven h 6 inches high, k 2 inches wide, and I9 inches deep. 
It is fitted with automatic temperature controls by means of which any 
desired thermal cycle could be imposed on the model. The furnace was 
equipped with a built-in loading frame suitable for the application of 
all basic types of loading. 

A special loading jig built for the investigation is shown in fig- 
ure 7« It consisted of a smooth circular shaft about I/2 inch in diam- 
eter passing through a pair of smooth, lubricated guide holes carefully 
alined so that the axis of the shaft was perpendicular to the base. The 
load was applied to the top of the shaft through a hard steel ball. This 
jig was found to give almost perfect vertical loading and the friction 
was negligible. 

A special jig was also built for oblique incidence. The frame of 
the jig could be rotated about a vertical axis through any desired angle 
which could easily be measured to one-tenth of a degree. The slice was 
mounted in the frame of the jig and the whole unit was placed in an 
immersion tank with a suitable mixture of Halowax and mineral oil. 



NACA TN 2822 


17 


The remaining equipment was standard apparatus in photoelastic 
laboratories the descriptions of. which are available in the literature. 


Test Procedure 

Model . - The sphere was machined from a cylinder of Fosterite which 
was previously annealed to remove initial stresses. The machining was 
carried out in a lathe. The rough cylinder of Fosterite was first turned 
to a true cylinder . In order to form the sphere a tool bit was set in 
a special Jig riding on the carriage of the lathe. This bit could be 
rotated about a vertical axis lying in the plane of the lathe centers. 

The cutting was performed by swinging the tool bit by heind around its ' 
vertical axis while the cylinder was rotating, and the radius of the 
sphere was slowly reduced by bringing the tool bit gradually closer to 
its axis. In this way it was possible to shape the complete sphere 
except for a relatively small nipple near the chuck. The final diameter 
of the sphere was 3.313 ± 0.002 inches. 

Loading .- The sphere was placed in the loading Jig and carefully 
alined for diametral compressive loading. The model was then heated to 
162° F in the electric furnace, the rate of heating being about 10° F 
per hour. A load of 172 pounds was then applied to the model. After a 
soaking period of about 2 hours the temperature was lowered at the rate 
of 4° F per hour to room temperature. The final diameter of the equa- 
torial plane was found to be 3.334 inches and the load axis was measured 
as 3.102 inches. Although relatively large deformations were developed 
in the loaded regions, the main body of the sphere was not notably dis- 
torted from its original shape. The stress pattern of the whole sphere 

in figure 8 shows that the loads and the stresses were rotationally 
symmetrical. 

. Slicing .- In preparation for slicing the center lines of all the 
slices were carefully scribed on the sphere using the flat spots in the 

loaded regions as datum planes. The slices were then sawed out roughly 

on a band saw to a thickness of about 3/8 inch. They were subsequently 
ground by hand to about l/8-inch thickness in most cases. Great care 
was taken to keep the slices symmetrical with respect to their center 
lines. 

Figure 9 shows the slicing plan. The first slice removed was 
parallel to the equator and midway between the equator nnd the load 
point. Then from the opposite side- of the sphere, a meridian slice was 
removed. Next a slice containing the equatorial plane was cut. Finally 
a slice parallel to the meridian slice and halfway out on the radius was 
removed. 



18 


NACA TN 2822 


Stress patterns and isoclinics .- The slices were moiaited in the 
oblique incidence jig and stress patterns at normal and oblique incidence 
were, recorded photographically. Typical stress patterns are shown in 
figures 10 to 15. Most of the normal incidence patterns show very few 
fringes. In order to obtain accurate data in these cases a Babinet- 
Soleil compensator was used to obtain the fringe-order distribution 
along the; lines of interest by point -by-point exploration. A black 
cross was made on the screen to mark the point on which attention should 
be centered. The model slice was then adjusted lontil the line of 
interest on the image ran true on the intersection of the cross as the 
straining frame was moved laterally. Then beginning at. the outer edge 
and moving the straining frame by a known amount after each observation 
the fringe order was obtained at a series of points along the line. 

From these data the curve of fringe-order distribution could be plotted. 

It may be noted that fractional fringe orders can also be obtained 
by the Tardy method of - compensation, the accuracy being comparable with 
that of the Babinet-Soleil compensator. 

Isoclinic lines were recorded by one of two methods . For the over- 
all picture the isoclinic .lines were photographed in most cases (figs. 16 
and 17). From the photographs averaged sketches were prepared and used 
in making the calculations. On several lines direct sketching of the 
isoclinics was used with attention being confined to the particular line 
of interest. Here the intersections of successive isoclinic lines with 
the line of interest were obtained at very short intervals, from which 
the distribution of the isoclinic parameters along the line could be 
plotted. This method was found to be accurate and considerably less 
time consuming than the photographic method. White light was used in 
all isbclinic work. 

In plane stress systems the isoclinic parameter at a point on a 
free boundary. is determined by the tangent to the boundary at the point. 
The isoclinic parameter thus changes from point to point along the 
boundary in general. This is not necessarily true for isoclinics of 
secondary principal stresses. In the slice parallel to the meridian the 
secondary principal stresses at the boundary consist solely of one normal 
stress Uz which is horizontal. The boundary is therefore a zero 

isoclinic and no other isoclinic may intersect the boundary at any point. 
The higher-order isoclinics therefore all lie within the boundary forming 
closed loops in this case, as shown in figure 17- 

Fringe value . - The fringe value of the material was obtained from a 

small cylinder about l/2 inch in diameter and 1^ inches long. * This was 

loaded in compression in the special loading jig used for the sphere and 
subjected to the same thermal cycle as the sphere. A portion of the 
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cylinder was machined away to leave a V shape as shown in figure 18. 

The resulting stress pattern was then photographed (fig. 19)* The 
V shape was used to make clearly visible the fringe of zero order 
occurring at the sharp edge of. the wedge. In the cylinder itself the 
first few fringes crowd together near the bo\indary of the cylinder and 
it is practically impossible to identify the zero fringe. From the 
stress pattern in figure 19 it was a simple matter to plot fringe order 
against position, which for the wedge described is a straight line 

(fig. 20). 

In this way the fringe order at the point where the plane faces of 
the V intersect the curved boundary of the cylinder could be accurately 
determined. At the same point the thickness h could also be measured 
accurately. With the measured fringe order and thickness the fringe 
value 2f of the material was found to be 3-0 psi per fringe, per inch 
compression. 


Re suits 

Interior stresses . ~ The stress distribution was obtained along six 
lines passing through the interior of the sphere. The lines are indi- 
cated in figure 21 by the letters A-A, B-B, C-C, D-D, E-E, and F-F. Thus 
the six lines are defined as follows: 

A-A intersection of meridian slice with equatorial slice 

B-B intersection of meridian slice with slice parallel to equator 

C-C intersection of s.lice parallel to equator with slice parallel 

to meridian 

D-D load line 

E-E intersection of equatorial slice with slice parallel to 

meridisui 

F-F intersection of slice parallel to the meridian slice with a 

meridian plane perpendiculeu:* to the above two slices 

In this problem advantage was taken of the symmetry of the sphere to 
eliminate the need for two models as discussed in the theoretical part 
of the report. In particular because of the rotational symmetry of the 
stresses one meridian slice could be taken to represent all meridian 
slices . 

Typical caJ-culation . - In order to make clear the application of the 
method the complete calculations for line C-C will now be given. The 
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evaluation of the stresses along this line requires all the generality 
which would be encountered in a body devoid of symmetry. The basic 
data for the determination of the stresses on this line are obtained 
from the stress patterns and isoclinics at normal incidence of the two 
slices defining the line C-C and from the stress pattern and. isoclinics 
of one of the slices at oblique incidence with rotation about an axis 
perpendicular to C-C and lying, in the plane of the plate. 

Because of the symmetry of the stresses along C-C it is necessary 
only to deal with half the length of the line. This half length was 
divided into 10 equal subdivisions. The two necessary auxiliary lines 
were drawn parallel to it in each of the two orthogonal planes and 
spaced the length of one subdivision apart. 

The first step is to obtain the distribution of the shearing 
stresses along line C-C and along the four auxiliary lines. This 
requires the determination of the distribution of the secondsiry principal 
stress differences and of the isoclinics along these lines. Figure 22 
shows the curves of n‘" and 0”’ for the slice parallel to the meridian, 
and figure 23 shows n" and 0” . for the slice parallel to the equator. 
With these data the required shearing stresses are computed at each 
division point of C-C. Thus, following equations ( 12) 

n^„ = i n’" sin 20"' (l6a) 

Hxz = I n" sin 20" (l6b) 

It will be noted that for positive values of z the shear systems nyz 
and nxz are positive. 

As noted in the theoretical part of the report the integration 
requires the use of the difference between the shearing stresses at the 
center of each subdivision. These differences are obtained from the 
curves representing the distribution of the shearing stresses just 
found. Figure 2h shows the curves of the shear differences for the two 
slices . 

The next step is to obtain the value of the normal stress n^ at 

each division point by ^ integration of one of the differential equations 
of equilibrium. The integration is carried out by approximation using 
summations to replace the integrals. The appropriate equation for 
line C-C is similar to equation (lla), that is. 
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Wj = Nc-r^y^i-i ^xzi (17a) 

c c 


Choosing Ax = Ay = Az the ratios of these quantities are unity. The 
signs of the ratios depend on the choice of axes and the direction of 
integration. The integration begins at the boundary and proceeds inward. 
Consequently Az is negative. The shear differences have been formed 
in such a way that Ay is positive and Ax negative. The final form 
of the equation is therefore 

(°z)j = (^z)c ^ ^yz -Y^ ^xz (1Tb) 

The signs of the shear differences are as shown by the curves. 

In order to start the integration the value of ^°z)c required. 

This must be determined from the boundary conditions and boundary fringe 
orders. From the fact that the boundary is unloaded it is evident that 
the principal stress normal to the boundary is zero. Inspection of the 
meridian slice shows that, excepting the loaded regions, the boundary 
stresses in the meridian section are also zero. The fringe order at the 
boTindary of the slice parallel to the equator is O.58 fringe tension, 
and the direction of this stress makes an angle of 55 - 6 ° with line C-C. 
The boundary value of is thus found from the equations of stress 

transformation as follows: 


(nzj^ = 0.58 cos255.6 = 0.185 (18) 

The expression for (nz)j therefore takes the form 

J 3 

(nz)j = 0.185 + YL ^yz - YZ ^xz (iTc) 

The integration is easily carried out in tabular form as shown in 
table I. , 

Once nz has been found the values of n^ and ny are found from 
expressions similar to equation (l3c), that is. 
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- n" cos 2ff>" (I9a) 

ny = nz - n'"cos 2ff)"' (l9b) 

This computation is shown in table II. 

The last step is to determine the remaining shearing stress 
system n^^y. This was obtained in this case from oblique incidence 

on the slice parallel to the- meridian with rotation through 45° about 
the Y-axis. Figure 25 shows the fringe order and isoclinic distribu- 
tions along C-C for this case. With these data and the known values 
of Uyz previously determined the required shearing stress component 

is found from an expression similar to equation ( 15 c). Thus 

"xy = \ “ey sin S^^y - Uy^ (20) 


The results of this computation are shown in table III. This completes 
the solution for all six stress components along line C-C. In order 
to determine the stress components in pounds per square inch it is 
necessary only to multiply the stresses in fringes by the proper fringe 
value of the slice. 

Using methods similar to those just explained the stress components 
for all six lines have been determined. With the exception of line D-D 
integration began at the boundary and proceeded inward. For line D-D 
the starting point was taken as the center of the sphere and integration 
proceeded upward. The starting value of. ny for this line was taken 

to be that obtained from line A-A. The results of these computations 
are shown in figures 26 to 32. At the center of the sphere the stress 
components were found to be Cy = -2.59F/A and “ 0.45P/A. 

These values may be compared with the stresses at the center of a disk 
\inder diametral compression which are Oy = -I. 9 IP/A and = 0.64 p/a. 

Contact stresses at load points .- As noted previously the loads on 
the sphere produced considerable local deformation which resulted in 
flattened areas at the poles. Upon observation of the meridian slice 
it was found that fringes and isoclinics were unusually clear right to 
the loaded boundary. It was therefore possible to continue the integra- 
tion all the way to the loaded boundary along lines normal to the areas 
•of contact and th'ereby to obtain approximations to the contact stresses. 
The stress components at the ends of line D-D represent the contact 
stresses at the poles, that is, at the centers of the loaded areas. In 
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order to determine, at least roughly, the actual distribution of the 
contact stresses over the loaded areas two additional slices normal to 
the loaded zones were cut from the remaining material of the sphere. 

The intersections of these slices with the meridian slice then define 
two lines parallel to line D-D which extend from line B-B to the loaded 
boundary. Starting values were taken from the stress distribution on 
line B-B previously obtained, and integration proceeded to the loaded 
boundary in the usual fashion. In this way the contact stresses at two 
points at different distances from the center were obtained. With these 
three points the distribution of the normal stresses on the contact 
surface could be pretty well determined. The shearing stresses acting 
on the surface of contact were found from the values of n* and the 
isoclinics 0^ in a meridian section in the region of contact. The 
results are shown in figure 33; the directions of the shearing stresses 
being from the poles outward. 

Checks on accuracy . - Two types of checks are available in this 
problem, static checks and checks between stresses on different lines. 
Static checks were made from the stresses on lines A-A and B-B and from 
those acting on a diameter in the surfaces of contact. Since these 
stresses are rotationally symmetrical the resultant force acting on the 
equatorial plane and on the pl^e containing line B-B parallel to the 
equator as well as on the plane of contact can be determined by integra- 
tion. From the stresses on line A-A the resultant load on the equatorial 
plane was computed as YjG pounds, which is 2.3 percent higher than the 
applied load of 172 pounds. The stresses on line B-B gave a resultant 
of 168 pounds which is 2.3 percent low. Lastly the resultaint of the 
normal stresses on. the surface of contact was found to be 17O pounds, 
or 1.2 percent low. 

Points 0* and 0]^' (fig. 2l) are common points on different lines. 

The stresses at these points can be determined from each line and the 
results compared. The stresses at 0* were found by integrating along 
line B-B and also by integration along the path AGO*. From line B-B 
the stresses were foimd to be n^ = -U.08 and nx = n2 = O.^lO. From 

the path AGO* they were found to be ny = -4,06 and nx = nz = G.42, 

At Gq’ the stresses are found from lines C-C and F-F. From line C-C 
the stresses were computed as n^ = -1.12, n^ - -O.36, n^ = G.^3, and 

nxy = 0 .62. From line F-F they were ny = -I.l4, n^ = -0.39, = 0.40, 

and n^y = 0 . 66 . The agreement in these values is seen to he quite good. 

It is also possible to compute the stresses on lines C-C and E-E 
from the stresses on lines A-A and B-B. Figure 3^ shows the recteingular 
stress components on lines C-C and B-B. The necessary equations for 
transformation are similar to the familiar equations for inclined planes 
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in plane stress systems, that is, 


ax + (Jy 
2 



cos 29 + T^y 


sin 20 


(21a) 


T 


0 



sin 20 


COS 20 


(21b) 


The stresses computed by transformation were compared with the stresses 
independently determined on lines C-C and E-E by integration. The com- 
parative values are shown in figures 28, 29, and 3I. The general agree- 
ment is seen to te excellent. 

Check by the Lame-Maxvell equations .- Lines A-A and D-D are lines 
of symmetry for the sphere. For these special lines the stresses can 
be computed by the method outlined by Jess op (reference I7) using the 
extension of the Lame-Maxvell equations to three dimensions. This 
computation has been carried out using fringe orders from the meridian 
and equatorial slices and the 5^ isoclinic from the meridian slice. 

The comparative values are given in figure 26 for line A-A and figure 30 
for line D-D. Inspection of the figures shows that in both cases the 
agreement is close. 


Analysis and Discussion 

The primary objective of the project under discussion was to develop 
a general method for solving three-dimensional problems photoelastically. 
In the theoretical part of this report such a method is described. The 
experimental work shows that the proposed method is practical. 

It is too early to draw broad conclusions regarding the general 
accuracy of the new method. However, the excellence of the static checks 
and the consistency of the results, as shown by the close cross checks 
between the results from the various lines, seem to indicate possibil- 
ities of high accuracy. Unfortunately there is as yet no theoretical 
solution available for this particular problem to furnish conclusive 
checks and a measure of the errors. Nevertheless there is a, reasonable 
degree of certainty that the major stresses are free from significant 
error. 

It must be pointed out that the stresses as found here represent 
the solution for a material for which Poisson *s ratio is 0.48, whereas 
most structural materials have Poisson *s ratios of about 0.3. This is 
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an inherent limitation of three-dimensional photoelasticity. However, 
as noted previously, the theoretical solutions available to date indicate 
that Poisson's ratio has only a small influence on the major stresses 
although the effect on the minor stresses may be pronounced (refer- 
ences 19 and 20). 

Although no theoretical solution is available for the sphere 
Hertz's solution can be used to check the contact stresses determined 
photoelastically . According to Hertz's solution the contact pressure 
is representable by the ordinates to a hemisphere erected on the contact 
surface. Further the Hertz theory predicts that the maximum contact 
pressure should be 1.5 P/Ac. Reference to figure 33 shows that the 
experimentally determined value of the maximum pressure is 1.53 P/Ac^ 
which is 2 percent high. When it is considered that the path of integra- 
tion used in determining the maximum pressure led along an equatorial 
radius to the center and thence up the load axis to the surface this is 
indeed a remarkable check. The general distribution of the pressures 
is also seen to be substantially correct. The proposed method would 
therefore seem to hold considerable promise for the determination of 
contact stresses . 

According to the Hertz theory the two principal stresses in the 
plane of the contact surface at the pole, Oy^ and should each 

equal — cTy. For Poisson *s ratio of I/2 this reduces to 

~ indicates that an isotropic point exists at the pole, 

and therefore the shearing stress is zero at this point. This is borne 
out by the photoelastic results. 

Along the circular boundary of the surface of contact the stress 

* /l - 2v\ P 

consists of a pure shear of the amount ^ according to the 

Hertz theory. If Poisson*s ratio is 1/2 this expression vanishes. 

Hence the normal and shearing stresses at the boundary should be zero. 

The photoelastic results are seen to give this value. 

The problem treated in this report has complete rotational symmetry 
which simplifies the experimental technique by e-liminating the use of 
two models- No problem has as yet been solved which requires two models. 
The use of two models will no doubt introduce complications, but no 
insurmountable difficulties are anticipated. However, further work must 
be done to demonstrate the effectiveness of the proposed method with two 
models. 

There remains also the possibility of using sub slices from the 
main slice as discussed in the theoretical part. The model used in this 
investigation was not large enough to make this procedure feasible 
although some attempts were made. This possibility also needs to be 
further explored. 
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It will perhaps also he desirable to repeat the solution of the 
sphere with smaller loads in order to reduce the relatively large local 
deformations . 


SUMMARY OF RESULTS 


The results from this investigation to develop a general method for 
three-dimensional photoelastic stress analysis may be summarized as 
follows: 


1. The method of strain measurement after annealing cannot be used 
with the materials now available in this country. 

2. A general photoelastic method of obtaining six stress components 
at any point of an unsymmetrical body arbitrarily loaded has been 
developed. This method does not depend on Poisson's ratio, althou^ the 
results reflect the physical constants of the model. 

3 . The new method is applicable in the plastic range of the model. 

4. The method shows possibilities for the determination of contact 
stresses . 

5 . The stresses existing in a sphere subjected to diametral com- 
pression have been determined with considerable accuracy. 

6 . At the center of the sphere the stress components were found to 
be Oy = -2.59 P/a and = Og ~ 0.|»-5P/A, where P is the load on the 

sphere and A is the area of the equatorial plane of the sphere. These 
values may be compared with the stresses at the center of a disk xinder 
diametral con 5 )ression which are dy = -I. 9 IP/A and = 0.64p/A. 

7 . Further work is needed to determine the full potentialities of 
the method when two models are used. Further work is also needed to 
determine the practicability of sub slices. 


. Illinois Institue of Technology 

Chicago 16 , 111 ., September 15, 1951 



NACA TN 2822 


27 


REFERENCES 


1. Prigorovsky, N. I., and Preiss, A. K. : Issledovanye Napryazionnavo 

Sostoyaniya Na Prozrachnykh Obyemnykh Modelyakh v Putchke 
Parallelnykh Luchey Polarizovanavo Svyeta (An Investigation of the 
State of Stress in Transparent Three-Dimensional Models in Parallel 
Beams of Polarized Ught). Izvestiia Akad. Nauk S.S.S.R. , Otdelenie 
Tekh. Nauk, no. 5, May 19^4^9, pp. 686-700. (See also M. M. Frocht's 
review of this paper in Appl. Mech. Rev., vol. 3, no. 2, Feb. 1950, 
p. 38 .) 

2. Oppel, G. : Polarisationsoptische Untersuchimg raumlicher Spannungs- 

und Dehnungszustande. Forsch. Geb. Ing.-Wes., Bd. 7, Nr. 5, 
Sept./Oct. 1936 , pp. 240-248. 

3 . Frocht, Max Mark: Photoelasticity. Vol. II. John Wiley & Sons, Inc.’, 

1948. 

4. Hetenyi , M. : The Fundamentals of Three-Dimensional Photoelasticity. 

Jour. Appl. Mech., vol. 5, no. 4, Dec. 1938 , pp. A-l49 - A-I 55 . , 

5 . Drucker , Daniel C., and Mindlin, Raymond D. : Stress Analysis by 

Three-Dimensional Photoelastic Methods. Jour. Appl. Phys . , vol. 11, 
no. 11, Nov. 1940 , pp. 724-732. 

6 . Drucker, D. C.: Photoelastic Separation of Principal Stresses by 

Oblique Incidence. Jour. Appl. Mech., vol. 10 , no. 3, Sept. 1943, 
pp. A -156 - A- 160 . (See also discussion by M. M, Frocht, Jour. 

Appl. Mech., vol. 11, no. 2, June 1944, pp. A-125 - A-126.) 

7 . Frocht, Max Mark: Photoelasticity. Vol. I. John Wiley & Sons, Inc., 

1941. 

8 . Weller, R. : A New Method for Photoelasticity in Three Dimensions. 

V Jour. Appl. Phys., vol. 10 , no. 4, April 1939, p. 266 . 

9 . Weller, R. , and Bussey, J. K.: Photoelastic Analysis of Three.^ 

Dimensional Stress Systems Using Scattered Light. NACA TN 717 

TQOO t J I f 


Id. Menges, Hermann Josef : Die experimentelle Ermittlung rdumlicher 

Spannungszustande an durchsichtigen Modellen mit Hilfe des 
Tyndalleffektes . Z.f.a.M'.M., Bd. 20 , Heft 4, Aug. 1940, pp. 210-217 

11. Hiltscher, R. : Polarisationsoptische Untersuchung des raumlichen 

Spannungszustandes im konvergenten Licht. Forsch. Geb. Ing.-Wes. 

Bd. 9 , Nr. 2 , March/April 1938 , pp. 9 I-IO 3 . * 



28 


NACA TN 2822 


12. Kuske, A.: Vereinfachte Auswerteverfahren ravunlicher spannungsoptischer 

Versuche. Z.V.D.I., Bd. 86, Nr. 35/36, Sept. 5, 19^2, pp. 54i_5l4.4. 

13* Leven, M. M., and Frocht, Max M.: Photoelastic Determination of 

Stresses on the Surface of Poppet Valves. Sci. Paper 158?} East 
Pittsburgh Res. Lab., Westinghouse Electric Corp., 1951- 

14. Leven, M. M.: Stresses in Keyways by Photoelastic Methods and 

Comparison with Numerical Solution. Proc. Soc. Exp. Stress Analysis, 
vol. VII, no. 2, 1949, pp. lUl-15l<-. 

15. Hetenyl, M. : A Photoelastic Study of Bolt and Nut Fastenings. Jour. 

Appl. Mech., vol. 10, no. 2, June 19^3^ PP* A-93 - A-100. 

16. Frocht, Max M.: Studies in Three-Dimensional Photoelasticity - 

Torsional Stresses by Oblique Incidence. Jour. Appl. Mech., vol. 11, 
no. 4, Dec. IShk, pp. A-229 - A- 23^. 

17* Jessop, H. T.: The Determination of the Separate Stresses in Three- 

Dimensional Stress Investigations by the Frozen Stress Method. 

Jour. Sci. Instr., vol. 26, no. 1, Jan. ,19^9, PP* 27-31* 

18. Kuske, A.: Das Kunstharz Phenolformaldehyd in,der Spannungsoptik. 

Forsch. Geb. Ing.-Wes., Bd. 9> Nr. 3> May/June 1938, pp. 139-1^9* 

19. Sternberg, E., and Sadowsky, M. A.: Three-Dimensional Solution for 

the Stress Concentration Around a Circular Hole in a Plate of 
Arbitrary Thickness. Jour. Appl. Mech., vol. I6, no. 1, March 19^9, 
PP* 27-38. 

20. Sadowsky, M. A., and Sternberg, E.: Stress Concentration Around a 

Triaxial Ellipsoidal Cavity. Jour. Appl. Mech., vol. I6, no. 2, 

June 1949, pp. 149-157. 



NACA TN 2822 


29 


TABLE I 

NUMERICAL INTEGRATION OF DIFFERENTIAL 
EQUATION (17c) FOR LINE C-C 


z 

Mean 

Mean 

nz 

COi' 

^yz 


0 



0.431 


0.012 

-0.009 


0.1 



0.410 


0.025 

-0 .024 


0.2 



0.361 


0.030 

-0.034 


Q.3 



0.297 


0.028 

-0.038 


O.h 



0.231 


0.021 

-0.038 


0.5 



0.172 


6.009 

-0.033 


0 

ON 



0.130 


0.001 

^ -0.024 


0.7 



0.105 


-0.004 

-0.0l4 


0.8 



0.095 


-0.019 

-0.001 


0.9 



0.113 


-0.060 

0.012 


1.0 



0.185 
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TABLE II 

CALCULATION OF NOR^ STRESS COMPONENTS FOR LINE C-C 




n" 

cos 29^" 

n'" 

cos 2 ( f )'" 

(1) 

nz - Dy 
( 2 ) 

“x 

( 3 ) 

(4) 

■ 


0.79 

1.000 

1.55 

1.000 

0.79 

1.55 

-0.36 

-1.12 

0.1 

o.4i 

0.79 

0.957 

1.52 

0.993 

0.76 

1.51 

-0.35 

-1.10 

0.2 

0.36 

0.78 

0.841 

1.43 

0.972 

0.66 

1.39 

-0.30 

-1.03 

0.3 

0.30 

0.77 

0.674 

1.30 

0.942 

0.52 

1.22 

-0.22 

-0.92 

0.4 

0.23 

0.76 

0.488 

1.16 

0.906 

0.37 

1.05 

BO 

-0.82 

I- ■ ■ 

0.5 

0.17 

0.74 


1.00 

0.863 

0.22 

■ 0.86 

-0.05 

-0.69 

0.6 

0.13 

0.70 


0.84 

0.827 

0.08 

0.69 

0.05 

-0.56 

0.7 

0.10 

0.66 

-0.035 

0.70 

0.809 

-0.02 

mm 

0:12 

-0.47 


0.10 

0.62 

-0.167 

0.56 

0.819 

-0.10 

0.46 

0.20 

-0.36 

QH 

0.11 

0.60 

-0.282 

0.4o 

0.866 

-0.17 

0.35 

0.28 

- 0.24 

1.0 

0.18 

0.59 

-0.375 

0.18 

' 

1.000 

-0.22 

0.18 

o.4o 

0 
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TABLE III 

ION OF SHEARING STRESS BY OBLIQUE INCIDENCE FOR LINE C-C 
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Figure 2 .- Necessary slices and directions of light for determination of 
surface stresses. L, direction of light. 
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Figure 5 «- Normal incidence and oblique incidence for different direc- 
tions of rotation. 




NACA TN 2822 


37 





Figure 6 .- Electric furnace and Control panel. 
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Figure 7.- Loading jig and model of sphere. 



Figure 8.- Stress pattern of whole sphere. 
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Figure 9-- Slicing plan for sphere. 1, slice parallel to equator j 
2, meridian slicej 3 ^ slice parallel to meridianj and 4, equa- 
torial slice. 
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Figure 10.- Stress pattern of meridian slice at normal Incidence. 




Figure 11.- Stress pattern of meridian slice for a rotation of 45 about 

Y-axis. 
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Figure I3.- Stress 



pattern of slice parallel to equator for a rota- 
tion of 45° about X-axis. 
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Figure 14.- Stress pattern of slice parallel to 

incidence. 



meridian at normal 



Figure 15- - Stress pattern of slice parallel to meridian for a rotation 

of about Y-axis. 
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Figure 16.- Typical isoclinic for meridian slice. 




Figure 17- - Typical isoclinic for slice parallel to meridian. 
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Figure 18.- Cross section of calibration member after machining of wedge. 



Figure 19.- Stress pattern of calibration member. 



Fringe Order 


NACA TN 2822 


45 



0 I 2 3 4 5 


Position , Arbitrary Units 

Figure 20.- Fringe order for calibration member , as a function of distance 

from edge of wedge. 



1^6 


NACA TN 2822 















Fringe Order 



Isoclinic Parameter, 4^ ,Deg 










Check from line B-B 


22 


53 



Figure 28.- Distribution of normal stresses along line C-C. 
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Fi<gure 30 .- Distribution of stresses along line D-D 
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Figure 31 - - Distribution of stresses along line E-E. 
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Figlire 3 ^-- Sketch relating stress components at a point on line B-B to 

to those on line C-C. 
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